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The formalism of the relativistic (or Dirac-) Brueckner approach in infinite nuclear matter is 
described. As nucleon-nucleon interaction the one-boson exchange potentials Bonn A,B,C and for 
comparison the Walecka model are used. The T-matrix is determined from the Thompson equation 
and is projected onto five covariant amplitudes. By the restriction to positive energy states an 
ambiguity arises in the relativistic Brueckner approach which is discussed here in terms of the 
pseudo-scalar and the pseudo-vector projection. The influence of the coupling of the nucleon via 
the T-matrix as an effective two-nucleon interaction to the nuclear medium is expressed by the 
self-energy. In particular we investigate the scalar and vector components of the self-energy for the 
different one-boson exchange potentials and discuss their density and momentnm dependence. We 
estimate the uncertainty of the self-energy due to the pseudo-scalar and the pseudo-vector choice. 
Usually the momentum dependence of the self-energy is thought to be weak, however, we find that 
this depends on the one-boson exchange potentials. For the Bonn potentials, in contrast to the auj- 
potential, the momentum dependence is strikingly strong above as well as below the Fermi surface. 
We compare with the results of other groups and study the effects on the equation of state and the 
nucleon optical potential. 
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I. INTRODUCTION 

In the relativistic (or Dirac-) Brueckner model the dynamical two-nucleon correlations in infinite nuclear matter 
are studied. The formalism is based on an effective quantum field theory for mesons and nucleons Q and thus ignores 
the underlying quark nature of the nucleon. In the Brueckner model the T-matrix (or Brueckner G-matrix) serves as 
an effective in-medium two-body interaction. It is determined by a self-consistent summation of the ladder diagrams 
in a quasipotential approximation (Thompson equation) to the Bethe-Salpeter equation, see, e.g., Ref. [||. The bare 
nucleon-nucleon (ladder) interaction is described by one-boson exchange potentials. In the present work we apply 
the Bonn potentials A,B,C of Refs. and for comparison also a simple uw-potential which is well known from the 
Walecka model Q . The nucleon inside the nuclear medium is treated as a dressed Dirac particle where the influence 
of the coupling to the surrounding nucleons via the T-matrix serving as an effective interaction is expressed by the 
nucleon self-energy. The self-energy has a Lorentz structure with large scalar and four-vector components and modifies 
the single-particle spectrum and the spinor wave functions. 

The relativistic Brueckner model has been developed in the last decades by various authors. At the beginnig of 
the 80’s the relativistic problem was considered by the Brooklyn group in first order perturbation theory. In 
the following the covariant formalism to a self-consistent treatment of the Thompson equation was developed by 
Horowitz and Serot and is in detail outlined in Ref. Q|. However, in these works the nucleon-nucleon interactions 
was described within the framework of the crw-model. Calculations with realistic nucleon-nucleon interactions have 
been performed by Brockmann and Machleidt and later on by ter Haar and Malfliet 00- A more rigorous 
derivation of the Brueckner approach in the framework of relativistic Green’s function can be found in Ref. @] . In 
the present work we refer to the formalism described in Refs. 0,0. 

The main difference between the relativistic Brueckner approach of Refs. and previous works is the use of 

medium dependent spinors and the requirement of self-consistency for both the spinor wave functions and the single¬ 
particle spectrum. The Dirac spinors of the nucleon in the medium differ considerably from the free spinors due to 
the large scalar and vector self-energies. These introduce an additional density dependence into the spinor matrix 
elements of the one-boson exchange potentials and thus into the T-matrix. This dynamical effect is absent in non- 
relativistics. Another new aspect of the modern Brueckner approach is that the inertial frames where the T-matrix 
is determined and where it is used have to be distinguished carefully. Due to the neglection of the retardation in the 
meson propagators the one-boson exchange potentials and the correspondingly T-matrix serve as quasi-potentials in 
the two-particle center-of-mass (c.m.) frame. Therefore the relativistic Thompson equation is solved in this system 
which can be done by numerical standard techniques, as, e.g., in the non-relativistic case of Ref. [T^ . On the other 
hand single-particle quantities like the self-energy of the nucleon are calculated in the nuclear matter rest frame, where 
the single-particle distribution, the Fermi sea, is naturally defined. Thus a covariant representation of the T-matrix 
is required to connect both frames. Even though, e.g., the single particle potential can be calculated with some 
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additional approximations avoiding this representation, as done in Ref. |^, the covariant technique is unavoidable in 
order to calculate the self-energy components themselves. 

In the relativistic Brueckner approach the self-consistent spinors are restricted to positive energy states. A different 
approach in the full Dirac space which is also based on one-boson exchange potentials is used in Refs. 10 . But in 
a full treatment one should further include renormalisation effects like the vacuum fluctuation contributions to the 
self-energy, as discussed in the context of the relativistic Hartree approximation in Refs. [|]Q|, and solve the full Bethe- 
Salpeter equation self-consistently. However, then the one-boson exchange potentials would be inappropriate. But 
these are still outstanding problems. In the standard Brueckner model which we discuss in the present work one gets 
around the problem of renormalization. The masses and coupling constants in the model are physical renormalized 
quantities and as a consequence one neglects the effective negative-energy states. However, the effective positive- 
energy spinors still contain admixtures of free negative-energy states because they are a superposition of free positive- 
and negative-energy spinors. To Lorentz transform the T-matrix as we have mentioned above the positive energy 
on-shell T-matrix is projected onto five invariant amplitudes. By the restriction to positive energy states an ambiguity 
arises which is discussed in terms of the pseudo-scalar and the pseudo-vector representation of the T-matrix in 
the Dirac space. From analogy to the one-pion exchange coupling in the nucleon-nucleon potential the pseudo-vector 
choice is more natural and further it is known from meson theory analysis that this choice suppresses the coupling 
to the negative-energy states. However, the difference between both choices also states a measure for the uncertainty 
concerning the determination of the self-energy components due to the restriction on positive energy spinors. One 
issue of the present work is to estimate this uncertainty. 

Thus, the aim of this work is to determine the self-energy components, to estimate their uncertainty due to different 
covariant representations of the T-operator and different nucleon-nucleon potentials and to examine the sensitivity 
of some commonly used observables, i.e., the equation of state and the optical potential on the self-energies. We 
are investigating in detail the density and the momentum dependence of the scalar and the vector components and 
compare our results to those of other groups Thereby the momentum dependence can just be determined 

in a first approximation because it is neglected in the self-consistency treatment of the Brueckner model. Commonly 
this approximation is justified by the assumption of a weak momentum dependence, at least inside the Fermi sea, but 
there have been also indications that this assumption depends on the chosen one-boson exchange potential ||^. We 
show here for the first time the momentum dependence over its full range, i.e., below and above the Fermi surface, 
and we find it in the case of the Bonn A,B,C potentials of the order of several 100 MeV, i.e., to be strikingly strong. 
In difference to the most publications in the field we focus our discussion to the self-energy components themselves. 

It is meanwhile well known that the nuclear matter saturation properties are much better reproduced in the 
relativistic than in the non-relativistic Brueckner approach, i.e., the predictions for the saturation point of a variety 
of nucleon-nucleon potentials are located along the Coester band which meets - in the relativistic case - the empirical 
area |^. The same also holds, e.g., for the nucleon optical potential. However, as we also find in the present analysis 
these commonly studied observables are not really sensitive to the self-energy components since they mainly depend 
on a cancellation of the components. The situation changes in the application of the Brueckner theory beyond infinite 
nuclear matter. When applying relativistic Brueckner results to the description of finite nuclei, e.g., in the framework 
of an effective (density dependent) relativistic mean field theory the self-energies start to play a decisive role. This 
fact is, however, most pronounced in the investigation of heavy ion collisions. Different from resting nuclear matter 
here the scalar and vector components are out of balance which is due to their different Lorentz transformation 
character and diverse observables are in particular sensitive to the magnitude of the fields |^. Furthermore the 
momentum space configuration can be highly anisotropic and rather corresponds to two counterstreaming currents 
of nuclear matter which can be described by two Lorentz-elongated Fermi ellipsoids. Approximate determinations of 
relativistic Brueckner mean fields for such configurations have been performed in p^ | and encouraging applications 
to heavy ion collisions in . Thus the density and momentum dependence of the nuclear self-energy components as 
the quantity which most prominant characterizes the medium effects is a question of particular interest. 

This paper is organized in the following way: in the next section we review the relativistic Brueckner app roach. 
There are still some technical differences in the standard Brueckner model between Ref. 0 and Refs. ||lo|-0. Here 
we follow the approach of Horowitz and Serot ra] w hich allows the use of standard techniques for the solution of the 
Thompson equation known from earlier works |]13| , P2[ . However, we use the modern one-boson exchange potentials, 
treat the isospin channels and also the real and imaginary part of the T-matrix separately and project on both, 
the pseudo-scalar and the pseudo-vector representation. In the following section we discuss our numerical results. 
The calculations with the crw-potential are compared with results from Ref. Q|. We further discuss the density and 
momentum dependence of the scalar and vector self-energy, the influence of pseudo-scalar and pseudo-vector choice is 
analyzed and we compare to the results of Brockmann and Machleidt ||^ and the Groningen group @0- Further the 
equation of state and the optical potential are shown and their sensitivity to the self-energy components is analyzed. 
In the final section a summary is given. 
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II. RELATIVISTIC BRUECKNER APPROACH 


In the relativistic Brueckner approach the nucleon inside the nuclear medium may be viewed as a dressed particle 
in conseqeunce of its two-body interaction with the surrounding nucleons. This problem is stated as a coupled set of 
three non-linear integral equations 


T = V + i J VQGGT, 

( 1 ) 

G = G° + G°EG, 

( 2 ) 

E = -i [ {Tr[GT]-GT). 

JF 

(3) 


The T-matrix is calculated in the ladder approximation of the Bethe-Salpeter equation ([^) and the bare nucleon- 
nucleon interaction is described by one-boson exchange potentials V. The intermediate nucleons are given by the 
two-body propagator iGG which is, however, usually replaced by an effective propagator |Q, i.e., the Blankenbecler- 
Sugar propagator or the Thompson propagator. In the present work the Thompson choice is applied. The Pauli 
operator Q accounts for the influence of the medium by the Pauli principle and projects the intermediate scattering 
states outside the Fermi sea. The dressed one-body Green propagator G is calculated via the Dyson equation (||) 
from the free propagator G°. The influence of the coupling to the surrounding nucleons is expressed by the self¬ 
energy S of the nucleon. In the Brueckner formalism this self-energy is determined by summing up the direct and 
exchange interactions with all the nucleons inside the Fermi sea F, see Eq. (^). Here the T-matrix plays the role of an 
effective medium dependent two-body interaction. In the relativistic Brueckner approach one introduces additional 
approximations for solving the equations (Q-||) which are sketched in the following, for further reading see, e.g., Refs. 


A. Self-consistent spinor basis 

The most general form of the Lorentz-structure of the self-energy compatible with translational invariance, her- 
miticity, parity conservation, time-reversal invariance and rotational invariance (in the nuclear matter rest frame) 
is 

S(fc) = Es(fc) - = Es(fc) - 7 oEo(fc) + 7 • ki;^(fc) (4) 

with a scalar part Eg and a vector part E^ = (Eq, kEt,) . Due to the Dyson equation i) the full Green function has 
the formal solution 

G-i ^ Qo-i _ ^ _ jijn[E(A;)] ( 5 ) 

where we have introduced the effective mass and the kinetic momentum 

M*(k) =M + ReS,(k), k*'^ = k>^+ ReS'^(k). (6) 

Re and Im denote real and imaginary part since (above the Fermi surface) the self-energy in general will be complex. 
Here we adopt the quasiparticle approximation |l^ , i.e., the Im [E] will be neglected in Eq. (0). This means that the 
decay width of the dressed nucleon is set equal to zero, resulting in an infinite lifetime of this ’quasiparticle’. In the 
relativistic Brueckner approach the full Green function G is replaced by its positive energy part. Thus the self-energy 
in Eq. (i is determined by the part proportional to the Eermi sea E, the so-called Dirac Green function Gd @0 

Goik) = -f M*{k)]2m5{k*'^ - M*^{k))Q{k*°)e{kF - |k|) (7) 

which eliminates the divergent contributions from the negative energy sea. However, in the form of Eq. ( 0 ) Gd 
is restricted to the description of nuclear matter at rest. In order to achieve general convenience the Eermi sphere 

0(fcF —|k|) has to be replaced by an Fermi ellipsoid ©(Ej? —with the Fermi energy Ep = y^(l -|- RyYkp'^ + 
and the boost 4-velocity = ( 7 , uy) |Q. In Eq. (^) the momenta of the quasiparticles are put on mass-shell, which 

is expressed by the ^distribution and thus = E*(k) = -\/k *2 + Hence the self-energy E(A:) = E(|k|,fcF) 

depends only on the three-momentum |k| and, of course, on the respective density ov kp. 
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Furthermore, in the relativistic Brueckner model one applies a mean field approximation to S, Eq. (jj), which allows 
a simple formulation of the self-consistency problem. The explicit momentum dependence of the self-energy which 
enters via the term kS„ can be dealt with by introducing the reduced kinetic momentum k*^ = k*^/{I £„) and 

the reduced effective mass M* = M*/(l -I- S.„). Further one neglects the momentum dependence of the effective mass 
{M* or M*). Thus, the nucleons fulfill a quasi-free Dirac equation 


UA(fc) =0 

and using the normalization of Ref. [S the self-consistent positive-energy spinors of helicity A are defined as 


u\{k) = 


lE*{k) + M* 
2M* 


1 

2A|k| 


E*(k)+M- 


X\ 


( 8 ) 


(9) 


with xa being a Pauli spinor. The Dirac spinors depend on the effective mass and thus on the nuclear density. 
The matrix-elements of the T-matrix and the one-boson exchange potentials V, Eq. (|^), are calculated with these 
spinors. In contrast to non-relativistic Brueckner theory this fact introduces an additional density dependence into 
the interaction which is one major reason for the great success of the relativistic treatment. This density dependence 
is mediated by the additional parameter M* which is fixed at a reference point, usually at |k| = kp. M* is obtained 
as the solution of the non-linear equation {reference spectrum approximation) 


M* = M + Ss(/cF, M*) - M*Ey{kF, M*) 


( 10 ) 


which follows from the definition given above. Self-consistency is achieved by determining for a start value of M* the 
medium-dependent spinors of Eq. , next the T-matrix from Eq. (|^ and finally the self-energy and the new M 
This procedure is repeated until convergence is reached. 


B. Covariant T-matrix amplitudes 


The T-matrix, in Eq. (|^, is most easily determined in the two-particle center-of-mass (c.m.) system while the 
self-energy, Eq. (^), is calculated in the nuclear matter rest frame. Therefore we will determine in a second step 
covariant amplitudes of the T-matrix. The Thompson propagator (and similar the Blankenbecler sugar propagator) 
projects the intermediate nucleons onto positive energy states and restricts the exchanged energy transfer by 6{k^) to 
zero. Thus equation (|^) is reduced to a three dimensional integral equation of the Lippmann-Schwinger type, the so 
called Thompson equation 10 


'r(p,q,a:)|c.m. = E(p,q) + 


(^ 


V{p,k) 


Q{k,x) 

E*^{k) 2E*{(i) - 2E*{k)+ie 


T(k, q,x). 


( 11 ) 


Here the Thompson equation is given in the two-nucleon c.m.-frame where p = (pi — P2)/2 = (Pi — P2)/2 is the 
relative momentum of the final states and similar q, k are the relative momenta of the initial and intermediate states, 
respectively. The Thompson propagator and similar the Blankenbecler-Sugar propagator imply that the timelike 
component of the momentum transfer in V and T is set equal to zero which is a natural constraint in the c.m.- 
frame, however, not a covariant one. Hence we solve the Thompson equation in the c.m.-frame in contrast to, e.g., 
the approach of Ref. 0. The starting energy in Eq. & is already fixed: a/s* = 2if*(q). The Pauli operator Q 
explicitely depends on the chosen frame, i.e., on the boost 3-velocity u into the c.m.-frame. With x we denote the set 
of parameters x = {kp,M*, |u|} on which the T-matrix actually depends. 

We solve the Thompson equation & for the on-shell T-matri x ( |p| = |q|) in the c.m.-system and thereby apply 
standard techniques which are described in detail by Erkelenz [ p2[ . We construct the positive-energy helicity T- 
matrix elements from the | JMLS” >-scheme. In the on-shell case only five of the sixteen helicity matrix elements are 
independent which follows from general symmetries, see Ref. |^. After a partial wave projection onto the \JMLS >- 
states the integral reduces to a one-dimensional integral over the relative momentum |k| and Eq. (0) decouples into 
three subsystems of integral equations for the uncoupled spin singlet, the uncoupled spin triplet and the coupled 


^From now on we omit the tilde because in the following we normally deal with M*, k*^. 
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triplet states. The Pauli operator Q is replaced by an angle averaged Pauli operator Q 0. Since the Fermi sphere is 
deformed to a Fermi ellipsoid in the two-nucleon c.m.-frame Q is evaluated for such a Fermi-ellipsoid. 


_ r 0 |k| < fc_ 

(5(|k|, a;) = {'yE*{'k) - EF)/u'y\'k\ for k_ < |k| < k+ 

If |k| > 


( 12 ) 


with k_ = \/kp — u’^Ep, k^ = ^{uEp + kp) and u = |u|. We are solving the integral equations by the matrix 
inversion techniques of Haftel and Tabakin 113]. Real and imaginary parts of the T-matrix are calculated separately 
by the principal-value treatment given in RefT]^ . Due to the antisymmetry of the two-fermion states we can restore 
the total isospin I of the two-nucleon system (1=0,1) with the help of the selection rule: (— = — 1 . From the 
five independent on-shell amplitudes in the \JMLS >-representation the five independent partial wave amplitudes in 
the helicity representation (for 1=0,1 and real and imaginary part separately) are obtained by inversion of Eq. (3.32) 
and then of Eq. (3.28) of Ref. [^. The summation over the total angular momentum J (we have taken J < 6 and 
for comparison J < 12) yields the full helicity matrix element 


E 


2J+1 

Att 


diy{e) < |p|A;A'|T^’^( a:)| |q|AiA 2 >=< pA^A'Il3|r(cr)|qAiA2ll3 > ■ 


(13) 


Here 0 is the scattering angle between q and p and A = Ai — A 2 ,A' = A] — A^. The reduced rotation matrices 
d\xi[d) are those defined by Rose [^. The matrix element (H) is actually independent of the third component of the 
isospin I 3 . The next step is the projection of the five independent helicity amplitudes onto five covariant amplitudes. 
Therefore the T-operator is expanded into the basis matrices in the Dirac spinor space 






(14) 


The subscripts (1), (2) denote the particle on which the matrix acts. This expansion, however, is not unique since 
we have specified the Dirac matrix structure of the T operator by its action only on positive energy states. The 
expansion of Eq. is given with g, pseudo-vector interaction ^Kys. It is defined as = 7 ^K^ with the 

transferred momentum at the vertex = p'f ~ Pi = P*/^ ~ P*^ where /, i indicates final and initial momentum. We 

will, however, also adopt a pseudo-scalar interaction where the fourth term in Eq. ( 0 ) is replaced by 
as it was also done by Horowitz and Serot Q. The pseudo-vector choice is, e.g., made by the Groningen group. 
This choice is supported by the argument that it agrees with the one-pion exchange coupling in the nucleon-nucleon 
potential and furthermore suppresses the coupling to the negative energy states which are neglected in the Brueckner 
approach @ 0 - In the present work we apply both variants since we want to examine the influence of the different 
choices on the self-energy. 

Next we take the helicity matrix elements of Eq. &■ On the r.h.s. the helicity states are acting directly on the 
matrix operators which are now abbreviated by Km for i=l,2 and Weis', V, T, P, A} 


< pA]A'H3|T(a:)|qAiA2H3 >= 0 < pA^A'|qAiA2 > T^’\\pl0,x). (15) 

m 

The helicity matrix elements on the r.h.s. are given explicitely, e.g., in Eqs. (2.11-2.17) of Ref. |2^. For on- 
shell scattering between positive energy states the matrix elements of the pseudo-scalar (PS) and the pseudo-vector 
(PV) matrix operators are identical This vertex equivalance follows immediately from the Dirac Eq. (H): 

^(p) ) 75 ^( 9 ) = u{p)^^u{q). Thus the covariant amplitudes are identical as well. The 

difference between both re pre sentations of the T-matrix in Dirac space enters only at the one-body level when 
calculating the self-energy 1^ . For fixed parameters (cc, |p| = |q|,0) equation 0) is a matrix relation between 
the five independent helicity amplitudes Tj (five of the sixteen amplitudes numbered with i = {A^A^; A 1 A 2 } = 1,.., 5) 
on the l.h.s. and the five unknown covariant amplitudes 

( 16 ) 

m 

This corresponds to Eq. (3.23) of Ref. ||^ with respect to our different normalization of the spinor basis, Eq. (^). The 
covariant amplitudes are determined by a matrix inversion of equation (li)- This, however, only has to be done for 
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two angles since for on-shell scattering in the two-particle c.m.-frame only two scattering angles appear, i.e., 0 = 0 for 
the direct and 6 = n for the exchange interaction. However, now the matrix becomes singular because two of the five 
helicity amplitudes vanish. Therefore it is useful to extract the leading angular dependence of the matrix Cim as well 
as of the helicity amplitudes Tj from the rotation matrices dj^yiO = 0, tt) in Eq. (O). Thus, the limit 6 * —> 0 , 7 r can 
be performed analytically. This method has been developed in appendix C of Ref. 0 and formula (^|) can be found 
there for the limiting angles 0 , 7 r (Eqs. (C.10,11)). Proceeding this way we determine the real and the imaginary 
part of the five covariant amplitudes for the direct and exchange contribution in both, the isospin singlet and triplet 
channels of the T-matrix = 0, tt). 


C. Self-energy components 


Now we are able to calculate the self-energy. Inserting the Green function, Eq. (|^, and the Dirac representation 
of the T-matrix, Eq. into the definition of the self-energy, Eq. (||), yields 


^a/3 — 


(fq 0(fcF-|q|) 
(27r)3 4£;*(q) 


M*1 




a(3 


ArpV _ rpS 


2 TV. 


OL0 


2T^- 


AT% - r| - 4T^ - 12rJ -f AT^ - ^ 




AM*"^ 


T 


4M*2 




- {A* - 


2q;ik*f^ - q*f^) 

AM*^ 


X 


tP 


(17) 


where the a, P denote both, Dirac (spin) and isospin indices. The direct and exchange parts of the T-matrix denoted by 
D and X are calculated from isospin sums of the isospin projection operators which yields -|- and 

Tjp = _X'™i=o -I- . The self-energy as a one-body operator in spinor space posseses components proportional 

to the unity matrix 1 and to 7 ^ as postulated in equation (^. From this relation the self-energy components Es,o,t; 
follow by taking traces 


Es{\k\;M*,kF)= ^Trp] 
Eo(|k|;M^fc^) = -iTr[ 70 S] 
E„(|k|;M^fcF) = -iTr[7•kS] 


d^q 

(^ 

d^q 

d^q 

(^ 


- |q|)-^^Ts(k,q;a;), 
e{kF - |q|)T(o)(k,q;a;), 

©(^F - |q|)^j^^^^r„(k,q;a;) 


where we have used the abbreviations 


rs(k,q;x) = - 


?'(o)(k,q;a;) = 


r„(k,q;a;) = - 


4T^ - T| - AT^ - 12TJ 


-4r, 


D 


-AT, 


D 


AT 


X 


T| - 2T^ - 2T^ 


T| - 2T^ - 2T^ 


k*^q*^^ - M* 
2 M *2 

E*(k) M*2 - /c*g*F 
E*{q) 2M*2 

|k| - k*q*>^ 


Tt 
t V 


T 


X 




2M* 


T 


X 


(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 


In the pseudo-scalar case [Q all the cofactors of Tj are replaced just by -1 which yields T(o) = Ty. The self-energy 
of a nucleon is calculated in the nuclear matter rest frame by integrating the effective two-body interaction over the 
Fermi sea. The 3-momenta k, q of a nucleon pair are given in this frame and we define P* = k -|- q = k* -|- q*. On 
the other hand the covariant amplitudes are determined in the two-nucleon c.m. system and thus depend on the 
respective quantities in that frame 


T^(kq,kq;x) = T"(|pe,„,. |, 0 = 0,x) + 3T'"i=i(|pc.„,,|, 0 = 0,cr), 

T^(kq,qk;x) = -T^^=°{\pc.m.\,9 = Tr,x) + 3T'"^=^(|pc.m.|, 61 = Tr,x). 


(24) 

(25) 


From the invariant mass s* = (E*(k) -|- E*(q))^ — P*^, the relative momentum in the c.m.-frame follows by |pc.m.| = 
y/ s*/4 — M*2. As already mentioned above, x abbreviates the set of parameters x = {kF, M* , |u|}, where the boost 
3-velocity (from the nuclear matter into the c.m.-frame) is now given by u = P*/vs* + P^. In the two-particle 
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c.m.-frame there exist only two possible scattering angles, i.e., 0 = 0 for the direct and 0 = tt for the exchange 
amplitude. Numerically we solve the Thompson equation © and determine the covariant amplitudes for a 
suitable two-dimensional range of the quantities Pc.m.,|u|. Using the azimuthal symmetry we calculate the self¬ 
energy, Eqs. (|l8[]20|) , by a two-dimensional integral over the Fermi sphere. Therefore for each nucleon pair we have 
to interpolate the covariant amplitudes at the respective values of Pc.m. and |u|. The resulting self-energy, taken at 
the Fermi surface |k| = kp, yields the new value of the effective mass M*, Eq. This value serves as the input 

for the next iteration. This procedure is repeated until a convergence of the effective mass is achieved. 


III. NUCLEAR MATTER RESULTS 

For the solution of the Thompson equation ([ll| ) we apply the Bonn A,B,C potentials of Refs. as the nucleon- 
nucleon interaction V. The potentials are calculated with the OBNNS code of R. Machleidt [^. The Bonn potentials 
are based on the exchange of the six nonstrange bosons with masses below 1 GeV (tt, 77 , p, w, (5, cr). For the pion a 
derivative pseudovector coupling is applied. The three parametrizations A,B and C differ essential in the ttNN form- 
factor and as a consequence in the strength of the nuclear tensor force. Actually we find only minor differences of a few 
MeV in the resulting self-energies, see below, and thus we present mainly the Bonn C results. For comparison with the 
results of Horowitz and Serot we also employ the simple uw-meson exchange of Ref. [Q with gl = 91.64, gl = 136.2 
and use the form-factor given therein. 


A. self-energy components 

1. density dependence 

First we want to demonstrate the reliability of our approach by comparison with the results of Ref. ||^. In Fig. 1 we 
show the self-consistent effective mass M *, Eq. as a function of the Fermi momentum kp which is a measure for 
the density g = 2/(37r^)fc^. Our results with the uw-meson exchange potential and the pseudo-scalar decomposition 
of the T-matrix reproduce quite well the values taken from Ref. [Q . The slight deviation for densities above saturation 
density (here kp = 1.42/m“^) is partly due to the different effective two-nucleon propagators. As mentioned above, 
we are using the Thompson propagator, Eq. (|^, while in Ref. @] the Blankenbecler-Sugar propagator is applied. 
Furthermore, in Ref. an additional approximation has been made. As discussed above, we are solving the Thompson 
equation for a two-dimensional grid of the relative momentum |p| and the c.m.-momentum |P| respectively the boost- 
velocity |u| and use this T-matrix for the determination of E. In Ref. 0 the |P|-dependence is replaced by an ’’average 
|P| approximation” for each |p| which causes deviations mainly at high density systems. For the Bonn C potential 
the effective mass is significantly different from that of the CTw-model. Here we applied the pseudo-vector choice. 

In Fig. 2 we show the density dependence of the self-energy, i.e., the scalar part Eg as well as the time-like —Eg 
and the space-like component —kpT,y of the vector self-energy. The self-energies are determined for the three Bonn 
potentials A,B,C and again we the pseudo-vector choice is applied. For the sake of a better comparison the sign of 
the vector self-energies is inversed and the dimensionless space-like component Ej, is multiplied hy kp. The scalar 
and also the vector self-energy component Eg,—Eg differ only slightly for the three Bonn potentials. The difference 
amounts about 5 MeV at saturation density (here kp = 1.35fm~^) and is still less than 10 MeV at kp = 1.90/m“^. 
In absolute magnitude Bonn C yields the weakest and Bonn A the strongest self-energies. The space-like vector 
self-energy —kpY,y is equal within line-width for the three Bonn potentials and is small in comparison to Eg, Eg, e.g., 
about 7% of Eg at saturation density. For further analysis we therefore concentrate on the large components Eg, Eg 
calculated with the Bonn C potential. 


2. pseudo-scalar versus pseudo-vector choice 

The two different representations of the T-matrix, i.e., the pseudo-scalar and the pseudo-vector case, represent 
an inherent ambiguity. This is, however, a common feature of all relativistic Brueckner calculations which specify 
the Dirac matrix structure of the T-matrix by using only positive energy T-matrix elements. In order to clarify the 
influence of this ambiguity in Figs. 3. and 4. we show the self-energy components Eg^g calculated exemplaryly with 
the Bonn C potential for both choices. In the pseudo-scalar case, see Fig. 3, the self-energies are larger by about 40-80 
MeV (decreasing with increasing density) than in the pseudo-vector case (Fig. 4). This difference may be interpreted 
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as a measure of the general uncertainty of the self-energy components arising from any particular projection of the T- 
operator on positive energy matrix elements only. This uncertainty is larger at low densities. To get more insight into 
the origin of this difference we decompose the self-energy in the two contributions stemming from the direct and the 
exchange T-amplitude respectively. E.g., in the pseudo-scalar case the direct part of the self-energy Sg is determined 
by the integral over Tq while the exchange part is given by the integral over the sum of all exchange-amplitudes 
-1/4T| - - 1 /4T|, see Eq. (^ ) and the subsequent comment. We find that in the pseudo-scalar case 

direct and exchange parts of Es,o are identical within less than 0 .5 MeV, see Fig. 3. However, in the pseudo-vector 
case the factor accompanying the 


TJ amplitude in Eqs. (|Tj^) reduces the exchange part by roughly 50%. For 
comparison in the crw-model, which lacks of pion exchange, the TJ amplitude is generally smaller and the effect, i.e., 
the reduction of the exchange part, amounts to roughly 25%. Of course also the direct part is different from that of 
the pseudo-scalar case since we consider the self-consistently determined fields. We conclude that the pseudo-vector 
choice in general reduces the total self-energy which is due to a considerable reduction of the exchange contributions 
while the direct contributions are even enhanced. 

In Fig. 5 we compare our results for Eg, —Eq obtained with the pseudo-vector choice and the Bonn C potential with 
the results of two other groups, i.e., the Groningen group 0 and the approach of Brockmann and Machleidt [|j. 
In Ref. the relativistic Thompson equation is solved by standard methods |0] directly in the nuclear matter rest 
frame instead of the two-particle c.m.-frame. Therefore the Thompson Eq. ( |ll|) is Galilei transformed to the rest 
frame with the c.m.-momentum P/2. The relativistic single particle potential U in the nuclear matter rest frame can 
be evaluated from this T-matrix by 


t/(|k|) = 


M* 

E*{k) 


Re< fc|E|fc>=^ 


M 


*2 


qsF 


E*{k)E*{q) 


Re < kq\T\kq — qk > 


(26) 


with the definition of E(|k|) as in Eq. (i- The advantage of this method is that one is not forced to Lorentz transform 
the matrix elements between the nuclear matter rest frame and the two-particle c.m.-frame and thus a projection 
onto covariant amplitudes is obsolete which simplifies the task considerably. However, now the self-energy components 
Eg, — Eq cannot be calculated directly as in the projection method of Eqs. (00) ■ Instead, one has to extract them 
from the single particle potential U. With the form of E as in Eq. (^), however, thereby neglecting E„ one finds 


u{m) = 


M* 


E*(k) 


Eg — En 


(27) 


i.e., the single particle potential results from a cancellation of the huge scalar and vector components. Using a mean 
field ansatz for Eg^Oj he-j neglecting their momentum dependence, one is able to determine them by a ht procedure to 
C/(|k|). In Ref. ^ the self-energy components are constructed by this potential-ht method. In Fig. 5 we compare these 
results with our calculations both for the Bonn G potential. In the density range considered the htted self-energies are 
generally smaller in absolute values by 50 up to more than 150 MeV, e.g., at saturation density they are reduced by 
about 75 MeV. In spite of this discrepancy both calculations have several features in common. E.g., the single particle 
potential, defined in Eq. (|^), is almost identical and as we will see later on the same holds for the equation of state. 
Concerning these two quantities the main difference between both treatments lies in the reference frame where the 
T-matrix is actually calculated while the projection method on covariant amplitudes is not essential. At moderate 
densities the T-matrix depends only weakly on the c.m.-momentum P and respectively on the chosen frame. Thus 
the numerical results show only minor differences. Another feature in common is that also the potential-fit method 
yields self-energies for the three types of Bonn potentials which are equal within a few MeV . 

In Fig. 5 we further compare to the results of the Groningen group [^. They also perform a full relativistic 
Brueckner calculation solving the Thompson equation in the two-particle c.m.-system and using a projection of the 
T-matrix onto a pseudo-vector representation. The difference to our method is on the one hand the application of 
a different one-boson exchange potential including different form-factors, the so called Groningen potential. On the 
other hand as - a more technical point - in Ref. |10(] the Thompson equation is solved in full momentum-spin space 
by the use of Fade approximants, for further details see, e.g.. Ref. |Q. For densities above saturation density the 
self-energies calculated with the Groningen potential are roughly the same as the ones obtained by us, however, in 
the low density range they are much smaller in magnitude. 


3. momentum dependence 


Up to now we have discussed the self-energy of a nucleon in nuclear matter obtained at the Fermi surface (|k| = kp), 
i.e., the density dependence of the self-energy. In the relativistic Brueckner theory the self-energy at the Fermi surface 
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and the T-matrix are determined self-consistently, as described above. The dressed nucleon propagators respectively 
the self-consistent spinor basis, Eq. (||) , entering into the calculation of the T-matrix are constructed with the effective 
mass, Eq. M* itself depends on the density dependent self-energies Es.o,t)(fcF)- Thus only the effective mass 

acts as a self-consistency or iteration parameter. In a next step the self-energies above and below the Fermi surface 
can be evaluated from Eqs. (p^pO|). This has to be considered just as the first iteration step in the determination 
of self-consistent momentum dependent functions Es_o,j;(|k|) and stands in contrast to the non-relativistic Brueckner 
treatment where self-consistency is required for the full positive single particle energy spectrum. However, these 
approximations are not based on the same footing. In the relativistic treatment we have an intrinsic momentum 
dependence due to the Dirac structure of the nucleon which introduces a momentum dependence in the single particle 
potential U even for constant values of see Eq. (|^). A further argument often given in connection with 

the reference spectrum approximation is that the momentum dependence of the self-energies in the relativistic 

Brueckner approach is found to be rather soft. However, in the present work we demonstrate that this is in general 
not the case but strongly depends on the special form of the one-boson exchange potential. In Fig. 6 we show 
the momentum dependence of the nucleon self-energies Ec g at saturation density calculated with the uw-exchange 
potential. Our results nicely agree with the results of Ref. ^ within 10 MeV which might be due to different technical 
details already mentioned in the discussion of Fig. 1. We find that decreases in the range from |k| = 0 to kp by 14 
MeV and Eg by 24 MeV, i.e., below the Fermi surface the self-energies are constant within a few percent. 

In Fig. 7 we show the momentum dependence of the self-energies for calculations with the Bonn C potential (pseudo¬ 
vector choice) and the crw-potential (pseudo-scalar choice), both at kp = 0.265 GeV. More strictly this is the real 
part of the self-energy since for momenta greater than the Fermi momentum the self-energy is complex. The results 
obtained with the Bonn C potential show a strong momentum dependence: over the entire momentum range up 
to 0.8 GeV corresponding to a single particle energy of around 300 MeV the self-energies decrease by nearly 60% 
with respect to the central value. The most pronounced change occurs in the region around the Fermi surface while 
deep inside the Fermi sea and far above the Fermi momentum the momentum dependence nearly saturates. For the 
CTW-potential the self-energies are rather weakly momentum dependent: over the full interval they decrease by about 
20%. For comparison we also show two results obtained by the Groningen group. In Ref. 0 the self-energies above 
the Fermi surface are calculated in a similar way as in the present work, see discussion of Fig. 5. The results obtained 
by the Groningen group are in general softer k-dependent than those of the present work obtained with the 
Bonn potentials. The results of Ref. show a momentum dependence as strongly as ours. However, in this case a 
comparison is not fully appropriate since these calculations include an additional self-consistent treatment of the pion 
polarisation. 


B. equation of state and optical potential 


In Fig. 8 we show the equation of state, i.e., the energy per particle E/A as a function of the nuclear matter density 
g, calculated with Bonn A,B,C. In the relativistic Brueckner theory the energy per particle is defined in analogy to 
the non-relativistic Hartree-Fock method as the kinetic plus half of the potential energy 


1 _ 1 A/* 

E/A = - ^< MA(k:)|7 • k -f M -I- -E{k)\ux{'k) > — - M 

^ k,A 


(28) 


1 

9 


27r^ 


((l + E.(|k|))A* 


E°(|k|)) 


— (E.dkDM^-E^dkDF'^) 


- M 


(29) 


with the self-consistent spinors |uA(k) > given in Eq. (^. By this way we have determined E/A by integrating the 
self-energy components of Eqs. (p^p0[) over the Fermi sphere, see Eq. (p^. In F^. 8 we compare to the results of 
Ref. 1^ where E/A is calculated directly from the single particle potential, Eq. (p^, using the expression (p^. For 
moderate densities the results agree with an accuracy better 0.5 MeV and even for kp = 1.9 fm~^ the deviation is less 
than 3 MeV. Firstly, this result is an additional test for the accuracy of our numerical treatment, i.e., the projection 
of the T-matrix onto covariant amplitudes which are subsequently used to determine the self-energy components and 
the equation of state. Secondly, these results again demonstrate that the equation of state is quite insensitive to the 
self-energy components themselves. Actually, as can be seen clearly from the definition Eq. (p^, the equation of 
state depends directly on the single-particle potential, i.e., the difference of Eg and Eg, Eq. (p^), and relatively weak 
on the effective mass via the kinetic energy part. Nevertheless, the relativistic effects are responsible in contrast to 
non-relativistic calculations that the minima of the equations of state for the three different Bonn potentials reveal a 
Goester band which meets the empirically found saturation area of nuclear matter. 

In order to demonstrate the influence of the two different choices, i.e., pseudo-vector and pseudo-scalar on the 
equation of state in Fig. 9 we show the respective equation of state, again for Bonn C. Actually the expectation value 
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< u\Y,\u > of the full self-energy operator, Eq. (|r^, is identical for both choices [^. However, as can be seen from 
Fig. 9 the equation of state turns out to be rather sensitive on the different choices. This significant dependence is 
only due to the self-consistently iterated effective mass which differs in both approaches. This effect deepens the 
equation of state in the pseudo-scalar choice by around 1 to 2 MeV and shifts the minimum towards higher densities. 
Thus the saturation point is shifted away from the empirical point towards the results of non-relativistic Brueckner 
calculations as, e.g., published in Ref. (H. This observation additionally supports the pseudo-vector choice. 

In view of the strong momentum dependence of the self-energy components for the Bonn potentials one might 
expect that the results depend sensitive on the actual momentum reference value where the selfconsistent effective 
mass is determined. As a consequence a selfconsistent treatment of the full momentum dependence might appear to 
be indispensable. Therefore we have also examined the influence of a different choice of the selfconsistent effective 
mass. In the reference spectrum approximation, see Eq. ( 0 ) , M* is determined at the Fermi surface which is the 
conventional treatment in the relativistic Brueckner approach. However, this choice is by no means unique. Thus we 
also determine M* by its mean value in analogy to Eq. ® 

M* = M+ %{M*) - (30) 


The self-energies T,s^v are not simply taken at |k| = kp but replaced by their values averaged over the Fermi sea 

/S,,„(|k|)0(A:^- IkDMVE^d^k 




/ 0(/cF - \k\)M*/E*d^k 


(31) 


These averages are Lorentz scalars. The extension of the reference spectrum approximation, Eq. (g0|), is especially 
meaningful with respect to the treatment of anisotropic nuclear matter, as done in Ref. |2^. Here it serves as a test 
of the influence of the momentum dependence. Due to the strong momentum dependence observed, e.g., in Fig. 7 
the mean value differs considerably from Es(/cf), be., by about 57 MeV at saturation density (with Bonn C). 
Respectively for the different effective masses we obtain a mean value M* = 528 MeV which is smaller than the 
selfconsistent mass M*{kF) = 566 MeV. But if we use the alternatively defined effective mass, Eq. (p0|), as iteration 
parameter we find for the final iterated mean value M* = 547 MeV. Thus the mean value is increased, i.e., it is shifted 
towards the previous self-consistent result or, in other words, the selfconsistent value is rather stable. The comparison 
of the self-energy components Tisfi^kp) in both iterations, see Fig. 2, shows just slightly deviations, e.g., about 10 
MeV at saturation density. And even more, the equation of state, shown in Fig. 9, remains nearly unchanged. These 
results demonstrate that the strong momentum dependence has only minor influence on the density dependence of 
the self-energy components and that the equation of state remains unaffected. 

Another quantity which is more sensitive on the momentum dependence of is the Schroedinger equivalent 

optical potential which a nucleon feels inside the nuclear medium. The optical potential is given by 


C/opt(|kU°) = E,(|k|)-^fc'^E^(|k|) 


E2(|k|)-E2(|k|) 


2M 


(32) 


It contains counteracting linear and quadratic terms in the self-energy components. For nucleons with |k| > kp the 
momentum space is not completely blocked for scattering processes and thus above the Fermi surface the self-energies 
and the optical potential become complex. The single-particle energy fc°(|k|) = — E° -|-F*(|k|) is a monotonic function 
of |k| which allows by inversion to determine f7opt(fc°) as function of the single-particle energy We want to stress 
that neglecting the momentum dependence of the fields E^yg) and neglecting completely E„ the optical potential 
Vopt would simply be a linear function of kP. However, this is not the case here. Real and imaginary part of the 
optical potential calculated for Bonn C (pseudo-vector choice) are shown in Fig. 10. The real part agrees well with the 
Groningen results of Ref. |11| while the imaginary part of the present calculations is about twice as large. Further we 
compare to the empirical optical potential Hama I, table 2 of Ref. For the full range of single-particle energies 
up to 700 MeV, shown in Fig. 10, our results deviate from the empirical Re C/opt by less than 10 MeV. Up to the pion 
production threshold at around 300 MeV we observe a remarkable agreement of the imaginary part of Uopt with the 
empirical data. For single-particle energies above this region the results appear to be no more fully reliable since then 
meson-nucleon resonances, e.g. the A(1232) resonance, should be taken into account. 


IV. SUMMARY 


In this work in particular the self-energy of a nucleon inside the nuclear medium has been studied in the framework of 
the relativistic Brueckner approach. Therefore the dependence of the Lorentz scalar and vector self-energy components 
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on the nuclear density and the nucleon momentum have been examined. As bare nucleon-nucleon interaction various 
one-boson exchange potentials have been employed, i.e., the Bonn A, B, C potentials and for comparison the auj- 
exchange potential. For the latter we are able to nicely reproduce the density dependence of the self-energies of Ref. 
0. The self-energy components themselves deviate for the three different Bonn potentials over a wide density range 
only by a few MeV and thus we restricted our discussion mainly to the Bonn C potential. From the comparison with 
our results we find that the potential-fit method of Ref. @ underpredicts the magnitude of the self-energies by about 
50 up to more than 150 MeV. 

Further we analyzed the influence of different possible covariant representations of the T-matrix. The positive 
energy on-shell T-matrix is projected on five invariant amplitudes and by the restriction to positive energy states an 
ambiguity arises which is discussed in terms of the pseudo-scalar and the pseudo-vector representation. Although the 
T-amplitudes of both choices are equal on the one-body level, i.e., concerning the respective self-energy components 
a significant influence of these representations is observed. This is due to the fact that in the pseudo-scalar case 
direct and exchange contributions to the self-energies are equal while in the pseudo-vector choice the exchange part is 
suppressed. For the Bonn potentials this suppression lies between 40 - 80 MeV and is most prominant at low densities. 
This difference between both choices has to be interpreted as an inherent measure for an uncertainty of the approach 
due to the restriction on positive energy spinors. 

Commonly the momentum dependence of the self-energy components is supposed to be weak. This assumption is 
also reflected in the self-consistency treatment in the relativistic Brueckner model which is managed by just one density 
dependent parameter, i.e., the self-consistent effective mass. Therefore we have carefully examined this momentum 
dependence which, however, can be done within the Bruckner scheme just in a first approximation. As in Ref. [Q we 
find that taking the crw-exchange potential the self-energies are found to be constant below the Fermi surface within 
a few percent. However, this result does not hold for the more realistic Bonn potentials. Here we observe a strong 
momentum dependence which may be due to the strong pion exchange. The self-energies decrease from their central 
value, i.e., at fc = 0 to the Fermi surface by about one fourth and further to A: = 0.8 GeV by more than one half. This 
striking momentum dependence makes the treatment of the self-consistency in the model somewhat questionable. But 
a test with a different defined self-consistent effective mass which is determined by a mean value over the Fermi sea 
and not taken at kp as usually done shows that the self-energy components are rather stable at not too high densities. 
Thus the self-consistency treatment in the model still appears to be justified up to two or three times saturation 
density. 

In this context the question raises which observables are sensitive to the self-energy components and their momentum 
dependence. It is remarkable that the full self-energy and, respectively, the single particle potential show (absolutely) 
a much weaker momentum dependence. As it is well known from mean field theory the large scalar and vector self¬ 
energy components contribute to these quantities with opposite signs and thus cancel each other to hight extent. The 
resulting observables, e.g., the single particle potential are smaller by nearly one order of magnitude. Here we found 
that this holds not only at kp but also over the full momentum range considered. Thus the momentum dependence 
cancels to a large extent. Also the equation of state depends mainly on the difference of the scalar and the vector 
self-energy and thus the absolute values of the single component are of minor importance. This can be seen most 
clearly by the use of just density dependent self-energy components determined with the potential-fit method of 
Brockmann and Machleidt Q which deviate from our results (at kp) by about 75 MeV. The resulting binding energy, 
however, only deviates by about 1 MeV. Similarly the self-energy components obtained in the pseudo-scalar and the 
pseudo-vector choice differ considerable but the equations of state differ by less than one to two MeV. On the other 
hand, small variations of the self-energy due to the different three Bonn potentials produce a significant change in 
the equations of state. Summarizing, the equation of state is not sensitive to the self-energy components themselves 
but to their difference. A bit more sensitive is another quantity, i.e., the optical potential. It depends not only on 
the cancellation of the scalar and vector self-energy but also contains quadratic terms. We find a good agreement of 
our results with the experimental data of Ref. | |2^ for, both, the real and the imaginary part up to a single particle 
energy of several hundred MeV. 

As a resume, the self-energy components in the relativistic Brueckner model are still affected with some uncertainty. 
For realistic one-boson exchange potentials we find them to be strongly momentum dependent. However, for this 
attribute we have to be aware of the limits of the model. One major success of the relativistic approach was that 
the equation of state comes close to the empirical saturation point. The equation of state is rather insensitive to the 
self-energy components themselves. Similarly, the nucleon optical potential is only little more sensitive. However, a 
precise knowledge of the size of each self-energy component is necessary, e.g., in the much more complicate situation 
of heavy ion collisions. 
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FIGURE CAPTIONS 
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Fig. 3: 
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Fig. 7: 

Fig. 8: 

Fig. 9: 


Fig. 10: 


Self-consistent effective mass M* in units of the bare mass (M=939 MeV) as a function of the Fermi momentum 
kp. The relativistic Brueckner results obtained with the troj-exchange potential (full line) are compared to those 
of Ref. [Q (diamonds). The effective mass for the Bonn C potential (dashed line) differs significantly from that 
of the (Tw-model. 

Self-energy components for the pseudo-vector choice, calculated with various Bonn potentials. The dotted, 
dashed and full lines represent the Bonn A,B and C potential. The scalar self-energy Ss (lower half-plane) 
and the time- and space-like vector self-energies —T,o,—kFY,y (upper half-plane) are shown as functions oi kp. 
For the three Bonn potentials the small space-like vector self-energies —kpY,y are equal within line width. The 
remaining curves (dashed dotted) involve a different definition of the self-consistent effective mass which is given 
in the text later on. 


Self-energy components — Eq calculated with the pseudo-scalar choice and the Bonn C potential. The curves 
represent the total self-energy components (full lines), their direct (dir.) and their exchange (exc.) parts which 
are practically equal (dashed dotted lines). 

Self-energy components E^, — Eq calculated with the pseudo-vector choice and the Bonn C potential. The curves 
represent the total self-energy components (full lines), their direct (dir.) parts (dashed lines) and their exchange 
(exc.) parts (dotted lines). 


Self-energy components E^, — Eq calculated with the pseudo-vector choice and the Bonn C potential (full lines). 
They are compared with the results of Ref. ^ (dashed lines) and of Ref. | (dotted lines). 


Self-energy components E^,—Eg at saturation density (kp = 0.28 GeV) as function of the nucleon momentum 
|k| < kp. The self-energies are calculated with the crtu-exchange potential (pseudo-scalar choice) and compared 
with values (symbols) taken from Ref. |^. 

Self-energy components E^,—Eg as function of the momentum |k| of the nucleon for kp = 0.265 GeV. The 
self-energies are calculated with the Bonn C (pseudo-vector choice, solid line) and the crw-exchange potential 
(pseudo-scalar choice, dotted line) and compared with results of Ref. [|^ (dashed line) and of Ref. Q (kp = 0.27 
GeV, dashed dotted line). 


Equation of state, i.e., energy per particle E/A as a function of the nuclear matter density p, for the three Bonn 
potentials A,B and C (from below to top: dashed line, dotted line and solid line). The symbols represent the 
respective results of Ref. Q . 

Equation of state for the Bonn C potential calculated in the pseudo-vector (solid line) and the pseudo-scalar 
(dashed line) choice. The remaining curve (dotted) involves a different definition of the self-consistent effective 
mass given in the text. 


Real and imaginary part of the Schroedinger equivalent optical potential Gopt as function of the single-paricle 
energy — M. The results for the Bonn G potential (real part: solid line, imaginary part: dotted line) are 
compared to the results of Ref. 0] (real part: dashed line, imaginary part: dashed dotted line) and to the 
empirical optical potential taken from Ref. ||^ (real part: diamonds, imaginary part: circles). 
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